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ABSTRACT 

A  semi-empirical  computational  model 
describing  the  design,  operation,  and  time- 
dependent  electrical  output  of 
magnetocumulative  generators  (MCGs)  has  been 
presented  in  previous  publications  by  the 
authors.  This  model  has  been  shown  to 
successfully  predict  the  behavior  of  a  number  of 
helical  generators  developed  by  laboratories 
around  the  world  over  a  number  of  years.  This 
paper  is  an  initial  attempt  to  describe  the  physics 
involved  in  the  time-dependent  behavior  of  the 
electrical  resistance  during  the  operation  of  these 
generators.  Understanding  this  parameter  is 
important  since  it  appears  to  be  dependent  on  the 
physical  dimensions  of  the  generator,  the 
detonation  velocity,  and  the  plasma  properties  of 
shock  heated  air.  The  early-time  decrease  in 
resistance  appears  to  be  related  to  the  rate  at 
which  the  plasma  is  formed,  while  the  late-time 
increase  in  resistance  is  related  to  the  transport 
properties  of  the  plasma. 
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I.  INTRODUCTION 


1950s.  Since  that  time,  a  number  of  other 
countries  have  started  research  programs  and 
there  is,  at  present,  several  significant  research 
and  development  programs.  One  of  the  most 
significant  programs  in  the  United  States  is  that 
of  Texas  Tech  University  to  understand  the  basic 
physics  and  the  limitations  of  MCGs. 

A  semi-empirical  model  for  helical 
MCGs  is  presented  in  Ref.  1 .  This  model  is 
physics  based  and  was  developed  to  describe  the 
physical  processes  that  take  place  in  the  MCG. 

It  was  used  to  successfully  describe  the  operation 
of  a  number  of  different  helical  MCGs  operating 
under  a  variety  of  conditions.  It  is  semi- 
empirical  in  that  it  is  based  on  a  single  adjustable 
parameter. 

One  important  characteristic  of  the 
helical  MCG  is  its  time  dependent  electrical 
resistance.  An  empirical  model  was  developed 
to  describe  this  resistance  as  a  function  of  the 
generator’s  properties.  This  paper  is  an  initial 
attempt  to  describe  those  physical  processes  that 
affect  this  electrical  resistance. 

A  brief  discussion  of  the  mathematical 
model  for  the  helical  MCG  is  presented  in  the 
next  section.  This  is  followed  by  a  detailed 
discussion  of  its  resistance. 


A  magnetocumulative  generator  (MCG) 
is  a  device  that  converts  the  chemical  energy  of 
explosives  into  electrical  energy.  This  is 
accomplished  by  locating  an  explosively  loaded 
armature  (or  liner)  coaxially  with  a  solenoid 
called  the  stator.  The  stator  is  activated  by  a 
seed  current  from  an  external  power  source. 
When  peak  current  is  flowing  through  the  stator, 
i.e.,  peak  magnetic  field  within  the  stator,  the 
explosive  in  the  armature  is  detonated  and  the 
metallic  armature  surface  starts  to  expand.  As 
the  detonation  wave  proceeds  down  the  length  of 
the  armature,  the  coils  of  the  solenoid  are 
sequentially  shorted  by  the  expanding  armature. 
This  process  leads  to  a  time  dependent 
amplification  of  the  output  current  from  the 
solenoid  into  a  load. 

Magnetocumulative  generators  were 
first  developed  independently  in  the  United 
States  and  the  Former  Soviet  Union  in  the  early 


II.  MATHEMATICAL  MODEL  FOR 
HELICAL  MCGS 

The  circuit  equation  for  an  MCG  is: 

+  RI  =  0,  (1) 

dt 

where  L  and  R  are  the  time  dependent  generator 
inductance  and  resistance,  respectively,  I  is  the 
electric  current,  and  t  is  time.  The  detonation  is 
initiated  at  t  =  0.  It  is  assumed  that  the 
inductance  and  resistance  behave  exponentially: 

L=L0e-t/T  (2) 

R  =  R0etlT,  (3) 

where  L0  and  Rr,  are  the  inductance  and 
resistance  just  prior  to  the  detonation  of  the 
explosive  and  i  is  a  characteristic  time.  This 
characteristic  time  depends  on  the  type  of 
explosive  used  and  the  dimensions  of  the  device. 
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The  value  of  i  for  the  Los  Alamos  Mark  IX 
generator  was  determined  empirically  to  be  48 


(IS. 

The  solution  of  Eq.  (1),  using  Eqs.  (2) 
and  (3),  is: 


(  ( 


I  =  /0  exp 
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where  In  is  the  seed  current.  Equation  (4) 
represents  the  behavior  of  the  generator  from 
time  t  =  0  to  a  time  tn  at  which  the  detonation 
process  ends  and  the  current  peaks  and  starts  to 
decrease. 


III.  MCG  RESISTANCE  MODEL 


The  resistance  can  be  written  in  the 


form: 


R  = 


<jAb 


(5) 


where  /  is  the  length  of  the  solenoid,  o  is  the 
electrical  conductivity,  and  AR  is  the  annular  area 
between  the  outer  surface  of  the  armature  and  the 
inner  surface  of  the  stator.  As  the  detonation 
front  proceeds  down  the  armature  and  the 
armature  expands  to  make  contact  with  the 
solenoid,  the  length,  /,  decreases: 

dl  l 

~  =  —  (6) 
dt  t 

and 


/  =  l0e 


-tlr 


(7) 

where  la  is  the  initial  length  of  the  stator.  Thus, 


where 


R  =  R0e 


-tlr 


R,  = 


ctAd 


(8) 

(9) 


This  corresponds  to  Eq.  (3). 

It  was  determined  empirically  [1],  that 

91/ 

R0=-^,  (10) 

■'o 

where  la  is  the  length  of  the  armature.  The 
constant,  91,  has  the  units  of  electric  field,  V/m. 
The  value  of  R0  for  the  Mark  IX  generator  was 
calculated  to  be  3.15  x  10'4  O. 

Relating  Eqs.  (9)  and  (10),  it  is  found 
that  the  electrical  conductivity  of  the  air  in  the 
space  between  the  armature  and  the  stator  is: 


a  =  _ljy[o — j-nho/rn  (in 

91/A 

The  area,  AR,  is  represented  by: 

AR=x(rs2  -ra2),  (12) 

where  rs  is  the  radius  of  the  stator  and  ra  is  the 
radius  of  the  armature.  The  values  of  the  various 
parameters  in  Eqs.  (11)  and  (12)  for  the  Mark  IX 
generator  are  [2]:  rs  =  17.8  cm,  ra  =  8.65  cm,  l0  = 
1.118  m,  4  =  1.73  m,  and  I0  =  0.5  MA.  Using 
these  parameter  values,  the  electrical 
conductivity  is  calculated  to  be  o  =  4.67  x  104 
mho/m. 

The  electrical  conductivity  can  also  be 
written  as: 


O  = 


ne 


m  y 


(13) 


where  ne  is  the  free  electron  density,  e  is  the 
electron  charge,  me  is  the  electron  mass,  and  v  is 
the  electron  collision  frequency.  This  expression 
assumes  that  the  electron  mobility  is  given  by: 


e 

//  = - . 

mev 


Writing  the  electron  collision  frequency  as 


V  =  no 


r  \1/2 

'  3kT  ' 


V  me  J 


(14) 


(15) 


and  using  Eq.  (13),  the  initial  free  electron 
density  is  found  to  be: 

ne  =  Gn^OmekTf\  (16) 

e 

where  na  is  the  ambient  atmospheric  particle 
density,  ac  is  the  electron  collision  cross-section 
with  air  particles,  k  is  the  Boltzmann  constant, 
and  T  is  the  temperature.  In  writing  Eq.  (16),  it 
was  assumed  that  the  degree  of  ionization  is 
relatively  low  and  that  the  electrons  are  in 
thermal  equilibrium  with  the  ambient  gas  in  the 
MCG.  In  most  cases,  the  gas  is  air. 

As  a  numerical  example,  consider  the 
ambient  gas  to  be  air  at  sea  level  density  ( na  = 
2.69  x  1025  m"3)  and  a  temperature  of  288.76°K. 
For  a  conductivity  of  4.67  x  104  mho/m,  it  is 
found  using  Eqs.  (16)  and  an  electron  collision 
cross-section  of  9.4  x  10"20  m2,  that  the  electron 
density  is  4.82  x  1023  nf\  This  corresponds  to  a 
1.8%  degree  of  ionization. 

As  pointed  out  in  the  previous  section, 
the  resistance  decreases  exponentially  out  to  a 
time  t0.  This  time  is  represented  by  the  formula: 
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t0=~,  (17) 

where  \,j  is  the  detonation  velocity.  The 
detonation  velocity  depends  on  the  explosive 
used.  As  an  example,  the  explosive  used  in  the 
Mark  IX  was  PBX9501,  which  has  a  detonation 
velocity  of  8.83  km/s  [2].  Thus, 

t0=  1.96X10"4  s. 

For  times  greater  than  t0,  it  has  been 
found,  empirically,  that  Eq.  (4)  must  be  modified 
to  the  following  form: 


7  =  /0exp] 

1 

-O 

of  ► 
1 

1 _ 

—  i 

l 

L 

rJ 

for  t  >t0,  where  b  =  1.8  x  105  s'1.  The  inverse 


of  b  is  a  second  characteristic  time,  5.56  ps.  It 
should  be  noted  that  Eq.  (18)  has  a  peak  value  at 
a  time  slightly  greater  than  t0.  As  an  example, 
again  using  the  Mark  IX  generator,  the  peak 
value  of  Eq.  (18)  occurs  at  t~  212  ps,  which  is 
approximately  16  ps  after  the  end  of  the 
detonation,  t0. 

The  form  of  Eq.  (18)  indicates  that  after 
t  =  to,  the  resistance  has  the  form: 


R  =  R0  exp 


b(t-t0) 


t 

X 


(19) 


for  t  >  t0. 


For  the  Mark  IX  generator,  the  resistance  has 
decreased  from  its  initial  value  of  3.15  x 

10"4  Ohms  to  5.31x10  6  Q  at  t  =  t0.  At  t  = 
212  ps,  where  the  current  reaches  a  peak  value, 
the  resistance  has  increased  to  6.78  x  10"5  Ohms. 
The  resistance  increases  to  its  initial  value  at: 

t=-\-  =  2.22xl0~4  s.  (20) 

1-  — 

bx 


Thus,  the  resistance  has  increased  by  a  factor  of 
59.3  in  a  period  of  26  ps  following  the  end  of  the 
detonation  and  continues  to  increase  very 
rapidly. 

IV.  SUMMARY 

This  paper  presents  an  initial  effort  to 
describe  the  physical  processes  involved  in  the 
time  varying  resistance  of  an  MCG.  The 
variation  in  the  resistance  was  shown  to  occur 


over  three  separate  time  regimes.  The  first  time 
period  is  that  following  the  onset  of  the  seed 
current,  but  prior  to  the  detonation  of  the 
explosive  charge.  The  second  time  period  is  the 
time  it  takes  for  the  detonation  wave  to  travel  the 
length  of  the  explosive  loaded  armature.  The 
third,  and  final,  period  is  the  time  following  the 
end  of  the  detonation  to  the  time  at  which  the 
current  in  the  load  has  been  reduced  to  zero. 

During  the  first  period  of  time,  the 
resistance  reaches  a  constant  value  that  is 
determined  by  the  amplitude  of  the  seed  current, 
the  dimensions  of  the  generator,  and  the 
surrounding  medium  (because  the  factor  91  has 
the  units  of  electric  field  and  because  of  the 
collision  frequency).  In  most  cases  studied  thus 
far,  the  surrounding  medium  has  been  air. 

The  resistance  during  the  second  time 
period  decreases  rapidly  with  a  characteristic 
time  that  depends  on  the  dimensions  of  the 
generator  and  the  type  of  explosive  used.  The 
rate  of  decrease  does  not  appear  to  be  dependent 
on  the  surrounding  medium. 

The  final  period  of  time,  following  the 
end  of  the  detonation  process,  the  resistance 
rapidly  increases.  The  characteristic  time  for  the 
increase  of  the  resistance  is  5.56  ps  and  is 
independent  of  the  generator  and  explosive 
characteristics.  This  would  seem  to  indicate  that 
it  is  dependent  on  the  surrounding  medium. 

Future  work  on  MCGs  will  include  an 
effort  to  determine  the  nature  of  the  dependence 
on  the  medium  in  which  the  interaction  is  taking 
place.  This  work  should  include  MCG  operation 
in  a  vacuum. 
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